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ABSTRACT 

A  theoretical  and  experimental  study  is  made  of  the  phenomenon  in¬ 
dicated  in  the  title.  It  is  shown  that  for  inertial  waves,  just  as  for 
gravity  waves,  there  are  discontinuities  in  the  tangential  particle 
velocities  at  the  interface  in  a  real,  two-component  liquid,  a  fact 
implying  the  existence  of  a  vortex  sheet.  For  the  case  where  the  two 
liquids  completely  fill  the  cylinder,  other  results  are  obtained  that 
are  analogous  to  those  for  gravity  waves.  In  particular,  if  the  liquids 
are  nearly  of  the  same  density,  there  are  two  sets  of  frequencies,  one 
set  characterizing  oscillations  of  the  liquid  mass  as  a  whole,  the  other 
set  characterizing  very  low  frequency  oscillations  at  the  interface. 

For  the  case  in  which  the  two  liquids  have  markedly  different  densities, 
there  are  again  two  sets  of  frequencies,  one  set  characterizing  oscilla¬ 
tions  of  the  inner  liquid  as  though  the  outer  liquid  were  a  solid  mass, 
the  other  set  characterizing  oscillations  of  the  outer  liquid  as  though 
the  inner  liquid  were  absent.  For  the  general  case,  physical  interpre¬ 
tations  are  difficult;  hence,  a  table  of  frequencies  versus  composition 
is  given. 


I 


i 

3 

% 

t 

t 

i 

s 

1 

\ 


i 

i 


4 

i 

f 

§ 

1 

? 

I 

I 

\ 

3 

3 

X 

i 

4 

? 

1 


3 


Preceding  page  blank 


«*•  i.V.v3.  MfVwi  wm**.  ami  ik  i  JK&iC 


TABLE  OF  CONTENTS 


Page 


ABSTRACT . 3 

LIST  OF  SYMBOLS . 7 

INTRODUCTION  .  9 

ANALYSIS  .  9 

A.  Mathematical  .  9 

B.  Experimental  .  .  .  .  . . .  .  19 

ACKNOWLEDGMENT . 20 

REFERENCES . 21 

DISTRIBUTION  LIST . 22 


K 


a 

b 

c 

d 

i 

J 


LIST  OF  SYMBOLS 

inner  radius  of  free  rurface 
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1 

K  20i/S 

2  _ 

K  ((2j  +  1)  iri/2d)/A  +  (2fl/8)2 

3 

p  pressure 
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I .  INTRODUCTION 


1 

S 

( 


Lamb  notes  some  interesting  facts  about  gravity  wave  motion  in 
inviscid,  incompressible,  superposed  liquids  of  different  densities, 
e.g.,  the  diminution  of  the  speed  of  a  wave  from  the  value  it  would 
have  in  a  one-component  liquid,  the  long  natural  period  of  oscillation 
of  the  interface  of  two  liquids  of  nearly  equal  density,  the  mathe¬ 
matically  possible  but  experimentally  unobserved  ''breaking"  of  the 
waves  at  the  interface,  the  independent  oscillations  of  the  two  liquids 
if  the  densities  are  markedly  different,  and  the  presence  of  a  vortex 
sheet  at  the  interface.  This  note  investigates  mathematically  and 
experimentally  whether  that  class  of  internal  waves  now  commonly 
termed  inertial  waves  (Bjerkness  ,  and  Fultz  )  exhibit  such  phenomena 
when  confined  in  a  right  circular  cylinder. 


II.  ANALYSIS 


A.  Mathematical 


3.  The  Governing  Equations.  Consider  Figure  1.  Following 

3 

Stewartson  ,  we  initially  assume  that  the  inviscid,  incompressible 
liquids  are  spinning  uniformly  as  a  rigid  body,  that  02b2  and  02c2 
are  much  greater  than  dg,  and  that  >  p^.  Then  in  this  steady  state, 

we  have: 


0  =  -  V 


rv.£ 
h  2 


(r2  -  a2)  ,  a  <  r  <  b 


0  =  -  ^  (b2  -  a2)  -  j-  (r2  -  b2)J  ,  b  <  r  <  c  (2) 

where  P  and  P  are  the  "static"  pressures  in  the  uniformly  spinning 
10  20 

liquids,  and  where  the  coordinate  system  has  the  angular  velocity 

n  =  (0,0,  n). 


References  are  listed  on  ‘page  21. 
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After  a  momentary  disturbar.ee  of  the  container  that  generates  a 
"perturbation"  velocity  u  =  (u,v,w) ,  the  linearized  Euler  equations  in 
the  two  regions  are: 


3  u  +  20  x  u 
It  "1 


V  j.  -  02  (r2  -  a 2) 

K  2 


3  u  +  2fi  x  u  =  -  V 
3t  “2  -  -2 


P  P  2 

L  2  2 


(b2  -  a2)  -  0^_  (r2  -  b2) 
2 


where  P^,  u^  and  P^,  u^  are  the  two  pressures  and  perturbation  velocities 

in  the  two  regions  after  the  disturbance.  Subtracting  (1)  from  (3)  and 
(2)  from  (4) ,  we  have 

3_uk  +  2fixuk  =  -V/  Pk  ~  Pko\=  -  V  jk  ,  k  ^  1,2  (5) 


where  PJ  and  P*  may  properly  be  termed  perturbation  pressures. 

4 

Again  following  Stewartson  ,  we  seek  normal  mode  solutions  of 

equation  (5) .  Hence,  letting  all  time-dependent  terms  have  the  same 

st  5 

functional  e  dependence,  one  gets  (Miles  ,  1959)  : 


*-r r-1- w 


Spk  [1  +  (20/S)2] 


From  the  continuity  equation  V  •  u^  =  0,  one  gets: 

,2  p.  .  (i|  .  ,)2  p,  .  o  (7) 

We  solve  equation  (7)  by  the  usual  separation  of  variables  process, 
getting : 


"  'z&wrzmsgmgFZ . 


P{  (r.e.z)  *  l 
\ 


^  k\(^ */  *  &  y  B"\(-?  r/T(¥2j 
1  **  * 


■[' 


cos  h  -j  z  +  D^  sin  h 


*■] 


(8) 


where  the  J's  and  Y's  are  Bessel  functions  of  the  first  and  second  kinds, 
respectively,  where  the  X*s  and  m's  are  separation  constants,  and  where 
A^,  B^,  C^,  are  integration  constants. 


2.  The  Frequency  Equation.  After  the  disturbance,  the  container 
has  only  the  angular  velocity  (o,o,fl);  hence,  the  boundary  condition 
on  all  solid  surfaces  is  u  •  n  =  0,  where  n  is  an  outward  directed  unit 
normal.  On  the  end  face  ,  then,  we  have: 


3Pk 


aZ 


] 


=  0 


2  =  ±  d 


(9) 


from  which  there  follows  =  0,  and 


\  =  *2  «  C2j  ♦  1)  i  -  ^T,  j  =  0,1,2,.  .  . 


(10) 


Hence,  redefining  our  constants,  we  rewrite  (8)  as: 


i  'll  ^  [*i  \  tK3r)  *  Bk  \  M  «  *  « 


^0mk 


HZ 

2d 


(ID 


where 


K3  .  (2,  ♦  1)  2d 


+ 1)  ii  /i  +/s£ 


(12) 


On  the  lateral  surface  r  =  c,  the  condition  u  •  n  =  0  gives: 


12 


3r  Sr  39 


r  =  c 


TK  J» 
2  I  3  “ 


’  (K  c)  +  K  _o  J  (K  c)-| 
m2  3  2  c  2  3 


+  b'  r 
2L 


K  Y'  (K  c)  +  K  %  Ym  (K  01  (13) 

o  TO  o  n  **'  TTl  I 


where 


_  2fii 
"  S 


In  section  B  we  outline  how  we  search  for  the  liquid  eigenfrequencies 
experimentally  by  adjusting  certain  physical  parameters  until  the  gyro¬ 
scope  undamps.  Inferring  that  this  undamping  is  due  to  resonance  between 

one  of  the  liquid  eigenfrequencies  and  the  nutational  frequency  of  the 

4 

gyroscope  (see  Stewartson  ,  1957),  one  can  determine  that  particular 

eigenfrequency .  Though  the  validity  of  this  procedure  is  well  documented, 

it  nevertheless  affords  a  determination  of  only  those  modes  for  which 

m  =1,  tor  one  can  show  that  coupling  between  the  liqud  and  shell 
2 

motions  occurs  only  for  that  value  of  m  .  Hence,  the  form  of  equation 

2 


(13)  pertinent  for  this  study  is: 

A'  fK  J»  (K  c)  +  \  J  (K  c)1 
2  L  3  1  3  C  1  3  J 

+  B'  fK  Y'  (Kc)+\y  (K  c)l 
2L  3  1  3  C  1  3  J 


=  0 


On  the  inner  free  surface  (at  r  =  a) ,  the  kinematical  boundary 

condition  is  ^-  =  0,  where  F  =  r  -  a  -  n  (9,z,t^  0,  and  n  (9,^i)is  the 

free  surface  elevation.  Hence,  to  first  order,  the  kinematical  boundary 
condition  is 


13 


? g^gpe^^g ?r^-.  v.iM^.»;V^ 


u  -  3rfl 

1  atl 


=  0  = 


r=a+n 


u.  -  s>] 


(16) 


r=a+n 


The  dynamic  boundary  condition  on  the  inner  free  surface  is  P  =  0. 

In  terms  of  P ' ,  this  is : 

1 

P*”]  =  -  P  I  s  -  o  n2an.  Hence: 

1-1  r=a+n  1®J  r=a+n  *  1 


SP'  +  p  fl2au  "1  =0 

1  l  1-1  r=a+n 


(17) 


which,  in  terms  of  (11),  is: 


A'  rj  (K  a)  -  K  nc  J'  (K  a)  +  K  \  J  (K  a)'l”| 

i  I  Bj  3  I  \  3  “j  3  2  a  1  3  J  J 

+  B*  fY  (K  a)  -  K  k  Y*  (K  a)  ♦  K  “i  Ym  (K  a)  11  *  0  (18) 

1  L  i  3  1  *•  3  “j  3  2  a  l  3  J  J 


where 


K  = 


fl2a 


‘  52  [-(f)2] 


(19) 


At  the  interface  r  =  b,  continuity  of  the  pressure  and  normal 
velocity  yield  two  dynamic  boundary  conditions  .  From  the  continuity  of 


the  normal  velocity  we  have: 

u  1  =  u  "1 

IJ  ...  2J 


,  or: 


r=b+n 


r=b+n 
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im  6 


m  **1  i 
1 

m 


fK  rA»  j* 

L 3  1 1  mi 


(K  b)  ♦  B'  Y'  (K 


1  m 


,  CK3b>} 


+  K  _j_  /A*  J  (K  b)  +  B*  Y  (K  b)\  “I 

2  b  1  1  \  3  1  }  J 

ie 

■  r-[Ks  {a;  j;  (Kab)  *  B;  y;  (K3b)} 

*  k2  i  {a;  (Kab)  *  b;  y>  «,">}  ] 


Hence,  m  =  1. 
1 


(20) 


From  the  equality  of  pressures  at  the  interface,  we  have 

f  P  ^  =  f  P  ^  ,  or  f  P 1  +  P  ^  =  fP'  +  P  >  .  Hence ,  to 

'  i'r=b+n  '  2'r*b+n  '  1  i^'r^b+n  '  2  2®'r=b+n 

first  order  in  n,  we  have:  (P!-PA)  •  bTj^  p_  -  p-)  ,  from  which 

x  r=>b+T) 

there  follows:  TsP1  +  p  02bu  =  SP’  +  p  02bu  "1  .  Substituting  for 

till  22  Zj-u+n 


P*  P1,  u  ,  and  u  ,  we  have: 
12  1  2 


A'  p 

r  (Kb)  . 

K  b  (k 

J'  (Kb)  + 

K  J 

(K  b))1 

i  L 

1  3 

la1 

3  1  3 

/  1 

3  J  J 

+  B' 

R  (K  b) 

-  Kb- 

f K  Y*  (K  b) 

+  K 

Y  (K  b)}l 

1 

L  i  3 

la 

l  3  1  3 

1  3  J  J 

=  A' 

[”j  (Kb). 

■  Kb/ 

K  J'  (K  b) 

*  K  J 

(K  b)>l 

2 

L  1  3 

la  >- 

3  1  3 

S* 

1  3  J  J 

+  B' 
2 

[Y.  (Ksb) 

-Kb- 

ijT 

fK  Y*  (K  b) 
l  3  1  3 

+  K 

S* 

Y.  «,«}] 

(21) 
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5-3* 


Equations  (15),  (18),  (20),  and  (21)  are  a  set  of  four  homogeneous 

equations  for  the  determination  of  A',  B',  A',  and  B*.  For  convenience, 

112  2 

we  re-write  these  four  equations  in  the  form 


A'  *  0  +  B'  •  0  +  A’J,  (c)  +  B'Jt,  (c)  =  0 
1  1  2  3  2  4 


(22) 


k'l  (a)  +  B'£  (a)  +  A'  •  0  +  B'  -0=0 
11  12  2  2 


(23) 


A'p  l  (b)  +  B*p  l  (b)  -  A'p  l  (b)  -  B'p  &  (b)  »  0 
123  124  213  21  4 


(24) 


A '£  (b)  +  B'Jl  (b)  -  A 'l  (b)  -  B'£  (b)  =  0 
11  12  2  1  2  2 

where  the  definitions  of  the  several  fi. * s  are  obvious. 


(25) 


From  the  above  set  of  equations  one  C2n  express  P'  and  P^  in  terms 


of  A'.  Then,  from  v  =  u  •  i„  = 
1  1-1-6 


L~S~  3r  r  30  J 


v  a  u  *  X-  = 
2-2-6 


[>  *  (f)2] 


r2o  9P’  i  9P’  i 

I  —  ^2.  _  _  J  #  one  can  show  that 


P2& 


v  j  #  v  *1  .  Hence,  as  is  the  case  for  gravity  waves  in  superposed 

l_>r=b  2-*r=b 

liquids1,  one  infers  that  for  inertial  waves  in  a  real  liquid  there 
would  be  a  vortex  sheet  at  the  interface. 

Equations  (22)  -  (25)  are  a  homogeneous  set,  and  a  necessary  condi¬ 
tion  that  A’,  B',  A'  and  B'  be  non- zero  is  that  the  determinant  of  the 
1  1  2  2 

coefficients  be  zero.  Hence: 
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0 

0 

*(c) 

3 

yc) 

1(a) 

1 

ya) 

0 

0 

c. 

=  0 

P  1  (b) 

Pi  (b) 

-  p  i  (b) 

"  Pi  (b) 

2  3 

2  4 

1  3 

1  4V  1 

yb) 

*  (b) 

2 

-  yb) 

-  i  o>) 

2 

(26) 

which  is  the  frequency  equation. 

We  first  consider  a  special  case  of  (26)  .  Let  the  two  liquids 
completely  fill  the  cylinder.  Then  a  =  0  and  equation  (26)  reduces  to 


£  (c) 
3 


yC) 


P  !  (b)  -  p  i  (b)  -  p  SL  (b)  =  0  =-  (p  -  p  )  4  (b)A  (b)i,  (c) 


2  1  1  3  h 


yb)  -  yb) 


«.  (b)  .  i.(c) 


[V, 


(bH  (b)  -  pi.  (b)i  (b) 
4  2  3  2 


We  now  consider  two  special  cases  of  equation  (27) . 

Case  I:  p  -  p 
2  *  l 

Then,  from  (27),  we  have: 

K 

i  (c)  =  K  J»(K  c)  +  -2- J  (K  c)  =  0  (28 

3  '313  C  iv  3 

which  is  Stewartson's  inertial  wave  frequency  equation  for  a  one- 
component  liquid  completely  filling  a  cylinder  of  radius  c. 

Alternatively,  from  (27)  we  have: 

yb)yb)  -  yb)yb)  a  yy)  [yyy) 

r  K  - 

-  Y  (Kb)  K  J'(K  b)  +r^J  (Kb)  =0  (2S 

1  3  L3  1  3  b  1  3  J 

Using  the  asymptotic  expansions  for  the  Bessel  functions,  one  can 
show  that  (29)  is  satisfied  if  ^  is  very  large.  However,  a  large  K 


17 


implies  a  small  3,  a  result  establishing  the  existence  of  low  frequency 
waves  at  the  interface. 


Hence,  we've  shown  that,  if  p  =  p  ,  there  are  two  sets  of  inertial 

1  2 

wave  frequencies,  one  set  for  the  liquid  as  a  whole,  the  other  for  the 
oscillations  at.  the  interface.  A  similar  result  holds  for  gravity  waves. 

Case  II:  p  »  p 
2  l 

Then,  from  (27),  we  have: 

K 

2  (b)  =  K  J'(Kb)  +r^J  (Kb)  =  0  (30) 

3  3  13  13 

This  would  be  the  Stewartson  frequency  equation  for  the  inner  liquid  if 
the  outer  liquid  were  a  solid  mass. 

Alternatively,  from  (27),  we  have: 


(bU  (c)  -  2  (c)  2  (b)  -  0 
14  3  2 


Cw  *  vb)  {y;(K3b)  few  *£vv°] 
[y;cv°  w][w)  -  Kitb)  hw  w}] 


(31) 


This  is  the  Stewartson  frequency  equation  for  a  one  component 
liquid,  with  a  free  surface  at  r  =  b,  in  a  cylinder  of  radius  c. 


Hence,  we’ve  established  that,  if  p  »  p  ,  there  are  also  two  sets 

2  1 

of  inertial  wave  frequencies,  one  characterizing  the  inner  liquid  as 
though  the  outer  were  solid;  the  other  characterizing  the  outer  liquid 
as  though  the  inner  were  absent.  A  similar  result  holds  for  gravity 
waves . 


For  the  general  case,  where  a  #  0,  and  pi  and  p^  are  neither  nearly 
equal  nor  far  removed,  equation  (26)  does  not  easily  lend  itself  to  such 


ag  ^m^f^rsp^r  -  —  — 


I 


physical  interpretations.  Hence,  we  merely  give  a  table  of  theoretically 
and  experimentally  determined  frequencies  versus  composition  for  a  par¬ 
ticular  density  ratio*. 


I 


E 


Table  I. 


p 

_L 

c 

b 

a 

p2 

d 

d 

.8 

3.127 

2.708 

0 

.8 

3.127 

2.211 

0 

.8 

3.127 

1.564 

0 

.8 

3.127 

2.752 

.989 

.8 

3.127 

2.319 

.989 

.& 

3.127 

1.783 

.989 

S  (Theory)  S  (Experimental) 

in  In 


931 

-  .936 

927 

-  .934 

938 

-  .942 

926 

-  .924 

914 

-  .920 

931 

-  .926 

B.  Experimental 

TVo  immiscible  liquids  are  placed  in  a  cylindrical  cavity  inside  a 

gyroscope.  Details  of  the  gyroscope,  supporting  apparatus,  and  the 

experimental  procedure  for  determining  the  linuid  eigenfrequencies  by 

resonating  one  of  them  with  the  nutational  frequency  of  the  gyroscope, 

4 

are  profusely  documented  elsewhere  .  We  merely  note  here  that  the 
inertial  waves  can  create  an  asymmetrical  pressure  distribution  in  the 
cavity;  and  the  above-mentioned  state  of  resonance  merely  ensures  that 
the  asymmetrical  pressure  distribution  will  "keep  step"  with  the 
nutational  motion,  thus  causing  an  overturning  effect  that  results  in  a 
growth  of  the  nutational  component  of  the  motion  of  the  gyroscope. 

Since  the  effect  involves  integrals  like 


o 


(cos  8,  sin  e)  d8 


that  are  zero  unless  m  «  1,  only  those  modes  of  oscillation  for  which 


*See  Section  B. 


m  *  1  couple  with  the  shell  motion.  Hence,  it  is  only  such  modes  that 
we  can  find  with  the  gyroscope. 

Since  viscosity  undoubtedly  has  an  effect,  it  was  essential  that 
the  two  liquids  of  different  density  have  small,  equal,  kinematic 
viscosities.  Fortunately,  the  Dow-Coming  Company  has  just  developed 
a  light  oil  that  has  a  kinematic  viscosity  (one  centistoke)  equal  to 
that  of  water,  yet  with  a  specific  gravity  of  0.8.  All  of  the 
experiments  were  run  using  water  and  this  oil,  and  the  results  are 
given  in  Table  I. 
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